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I.  INTRODUCTION 


The  study  of  the  motion  and  stability  of  a  liquid-filled  projectile 
requires  a  knowledge  of  the  internal  motion  of  the  liquid  which  determines  the 
moment  exerted  by  the  liquid  on  the  projectile.  The  motion  of  the  spinning 
liquid  can  be  divided  into  two  categories:  (1)  solid  body  rotation,  the 
condition  attained  for  "large"  time,  and  (2)  spin-up,  the  transient  state 
before  solid  body  rotation  is  attained.  Investigations  have  been  directed 
along  two  lines:  (1)  the  determination  of  the  frequencies  and  decay  rates  of 
inertial  oscillations,  and  (2)  prediction  of  the  moment  on  a  spinning  liquid- 
filled  cylinder  executing  angular  motion.  These  are  referred  to  as  the 
"eigenvalue  problem"  and  "moment  problem,"  respectively. 

For  solid  body  rotation  rational  theories  have  been  developed  and  applied 
both  to  the  eigenvalue  problem1  and  to  the  moment  problem  with  restrictive 
conditions. 2 3 4* 3,4  The  works  of  Stewartson5  and  Wedemeyer6  form  the  basis  of 
these  efforts.  Clearly,  to  include  the  transient  spin-up  effects  requires  a 
non-trivial  extension  of  those  works.  Eigenfrequencies  and  decay  rates  have 
been  obtained  for  spin-up;7’8  the  solutions  did  not  satisfy  the  no-slip 
condition  at  the  endwalls.  In  Reference  7  it  was  shown  that  calculated 
frequencies  compared  with  experimental  results  to  within  a  few  percent,  but 
decay  rates  could  be  off  by  as  much  as  a  factor  of  two.  If  an  endwall  cor¬ 
rection  for  solid  body  rotation  is  applied  during  spin-up  (obviously  an  ad  hoc 
correction)  the  prediction  of  decay  rate  is  improved  for  late  time  but  not 


1.  C.  W.  Kitchens ,  Jr.,  N.  Gerber ,  and  R.  Sedney ,  "Oscillations  of  a  Liquid 

in  a  Rotating  Cylinder:  Part  I.  Solid  Body  Rotation  "  BRL  Technical 

Report  ARBRL-TR-02081 ,  June  1978.  (AD  A057759) 

2.  N.  Gerber ,  R.  Sedney ,  and  J.  M.  BartoSj  "Pressure  Moment  on  a  Liquid- 

Filled  Projectile :  Solid  Body  Rotation ,"  BRL  Technical  Report  ARBRL-TR- 
02422 ,  October  1982.  (AD  A120S67) 

3.  C.  H.  Murphy ,  "Angular  Motion  of  a  Spinning  Projectile  with  a  Viscous 

Liquid  Payload "  BRL  Memorandum  Report  ARBRL-MR-03194 ,  August  1982. 

(AD  A1 18676 )  (See  also  Journal  of  Guidance ,  Control ,  and  Dynamics ,  Vol.  6> 
July-Auoust  1983 ,  pp.  28C=GTS6~. ) 

4.  N.  Gerber  and  R.  Sedney ,  "Moment  on  a  Liquid-Filled  Spinning  and 

Nutating  Projectile:  Solxd  Body  Rotation "  BRL  Technical  Report  ARBRL- 
TR-02470 ,  February  1983.  (AD  A12S332) 

5.  K.  Stewart8on3  "On  the  Stability  of  a  Spinning  Top  Containing  Liquid  ," 

Journal  of  Fluid^  Mechanic*} ,  Vol.  5,  Part  4,  1959. 

6.  E.  H.  Wedemeyer ,  "Viscous  Corrections  to  Stewartson' s  Stability 
Criterion  ”  BRL  Report  No.  1287 ,  June  1966.  (AD  489687) 

7.  R.  Sedney  and  N.  Gerber ,  "Oscillations  of  a  Liquid  in  a  Rotating 

Cylinder:  Part  II.  Spin-Up "  BRL  Technical  Report  ARBRL-TR-02489 ,  May 
1983.  (AD  A129094) 

8.  C.  H.  Murphy ,  "Moment  Induced  by  Liquid  Payload  During  Spin-Up  Without  a 
Critical  Layer  "  BRL  Report  in  preparation.  Also  AIAA  22na  Aerospace 
Sciences  Meeting ,  Reno ,  Nevada ,  January  19843  AIAA  Paper  No.  84-0229. 
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early  time.  In  this  paper  we  calculate  the  contribution  of  the  pressure  to 
the  moment  during  spin-up  of  a  liquid  that  fills  a  right  circular  cylinder 
undergoing  coning  motion  at  constant  frequency  and  constant  yaw  angle.  The 
contribution  of  the  wall  shear  will  be  treated  in  a  later  report. 

Two  coordinate  systems  will  be  used.  The  first  is  an  inertial  system; 

i.e.,  the  unyawed  reference  frame  shown  in  Figure  1,  with  rectangular 

2 

coordinates  (x,  y,  z)  and  cylindrical  coordinates  (r,  e,  x),  where  r  =  (y  + 

2  1/2 

z  )  '  and  G  =  arctan  (z/y).  The  second  reference  frame  will  be  defined 

later. 

Here,  as  in  previous  work*2’3*4  we  assume  the  yaw  to  be  sufficiently 

small  so  that  a  linearized  analysis  is  applicable;  i.e.,  the  flow  may  be 

considered  the  sum  of  a  basic  unperturbed  flow  and  a  perturbation  flow.  Also 
we  specify  that  the  perturbed  flow  variables  have  the  form* 


et[(l  -  ie)  x  it  -  0]  fn  (r>  x)>  (1.i) 

where  t  is  time,  x  and  e  are  non-dimensional . nutational  frequency  and  yaw 
growth  rate.  The  x  is  non-dimensional i zed  by  the  spin  of  the  cylinder;  r 

and  x  are  non-dimensional i zed  by  a,  the  cross-sectional  radius  of  the 
cylinder.  This  form  of  solution  renders  the  problem  time-independent  (aside 
from  the  quasi-steady  spin-up  perturbation  assumption,  to  be  discussed  later); 
and  it  implies  that  there  are  no  transients  present  between  the  angular  motion 
of  the  cylinder  and  the  corresponding  response  of  the  flow.  One  would  not 
expect  this  assumption  to  be  valid  early  in  the  spin-up  history  of  the 
rotating  liquid,  but  at  least  it  should  be  applicable  when  solid  body  rotation 
has  been  attained.  Justification  for  this  assumption  will  depend  on 
comparison  of  the  results  with  those  of  measurements. 

The  second  coordinate  system  used  to  describe  the  projectile,  Figure  1, 

is  the  x,  y,  z  non-spinning  system  which  has  the  x  -axis  along  the  projectile 

axis  of  symmetry;  the  y  and  z  axes  are  omitted  for  clarity.  The  x  =  0  and 

x  =  0  values  are  located  at  the  midplanes  of  the  unyawed  and  yawed  cylinders, 

respectively.  The  x  -axis  is  nutating  about  the  x-axis  with  the  angle  K^(t); 

the  pivot  point  lies  at  the  midplane.  The  components  of  the  projection  in  the 

y,  z  plane  of  a  unit  vector  lying  on  the  x  -axis  are  denoted  by  ny£  and  n^, 

respectively.  It  is  convenient  to  combine  the  two  components  of  yaw  into  a 
single  complex  variable 


C  =  -(ny£  +  i  nZE). 


(1.2) 


The  nomenclature  here  is  that  prescribed  in  Reference  3  and  used  in  References 
2  and  4. 


*Defintion8  of  quantities  are  given  in  the  LIST  OF  SYMBOLS  Section. 
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The  projectile  motion  is  specified  as: 


i  f$t 
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(1.3) 


where 


Kj  =  KQeft^t>  4>j  =  T^t,  f  =  (1  -  ie)  t  . 


(1.4) 


Here  K0  is  the  magnitude  of  the  yaw  at  time  t  =  0;  the  yaw  grows  when  cx  >  0. 
Also  <|>j  is  the  angular  orientation  of  the  x-axis  in  the  x,  y,  z  system  as 
shown  in  Figure  1. 

The  specified  projectile  motion  given  in  Eq.  (1.3)  was  found  in  gyroscope 
experiments  to  be  applicable  to  a  portion  of  the  coning  motion  history  in 
which  the  fluid  was  practically  spun  up.2’4  However,  this  simple  motion  is 
not  expected  to  occur  very  early  in  flight  while  the  spin-up  process  is  domi¬ 
nating  the  flow;  thus,  our  model  will  not  simulate  the  free-flight  situation. 
We  turn,  instead,  to  forced  coning  motion  and  fix  t  and  e  as  input;  in  par¬ 
ticular,  we  choose  e  =  0,  implying  no  yaw  growth.  Emphasis  will  be  placed  on 
investigating  resonances  between  the  angular  motion  and  free  oscillations  of 
the  liquid. 

After  the  t,  o  portion  of  the  solution  in  Eq.  (1.1)  is  factored  out,  the 
problem  is  posed  in  the  r,  x  plane.  We  attempt,  as  in  References  2  and  4,  to 

solve  it  by  a  modal  analysis,  i.e.,  by  separation  of  variables.  For  solid- 
body  rotation,  we  succeeded  in  obtaining  a  rational  solution;2  this  solution 
required  a  corrected  endwall  boundary  condition  for  axial  velocity  in  order 
for  it  to  satisfy  the  no-slip  condition.  Because  of  the  simple  basic  flow 
(fluid  angular  velocity  about  the  axis  =  $>  everywhere)  the  corrected  boundary 
condition  is  a  simple  one,  containing  a  complex  constant,  5c,  proportional  to 

Re~*^,  where 


(1.5) 


is  the  Reynolds  number  and  v  is  the  kinematic  viscosity  of  the  liquid. 

In  the  present  problem  the  basic  flow  is  the  spin-up  flow,  for  which  a 
corrected  endwall  condition  is  not  available.  The  Ekman  layers  (see  Reference 
7)  should  be  included  in  the  basic  flow  but  are  not  in  the  present  analysis. 
The  effect  that  including  them  would  have  on  a  separation-of-variables  ap¬ 
proach  is  not  known  at  this  time. 

Murphy8  has  proposed  a  method  of  calculating  the  liquid  moment  which  has 
the  intent  of  bypassing  the  difficulty;  namely,  to  apply  an  endwall  boundary 
condition  on  the  axial  flow  identical  in  form  to  that  for  solid  body  rotation, 
with  the  difference  that  the  Re  in  the  expression  for  the  6c  previously  men¬ 
tioned  is  replaced  by  an  "effective  Reynolds  number,"  Re^.  This  Reynolds 

number  is  obtained  by  averaging  over  a  meridional  plane  some  dynamic  quantity 
of  the  basic  spin-up  flow  which  varies  with  r,  e.g.,  azimuthal  velocity  or 
angular  momentum.  In  this  model  the  "corrected"  endwall  boundary  condition 
for  axial  flow  is  fitted  in  a  least  squares  sense. 
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This  procedure  is  not  a  rational  approximation,  defined  as  one  in  which 
the  next  term  in  the  approximation  (or  the  error)  can  be  calculated  in 
principle  or  estimated  in  an  order  of  magnitude  sense.  Therefore,  it  must  be 
considered  an  ad  hoc  approximation.  The  solution  satisfies  an  incorrect 
boundary  condition  approximately;  there  is  no  way  of  estimating  the  error  in 
the  moments.  Obviously,  the  approximation  is  more  valid  for  small  departures 
from  solid  body  rotation.  This  ad  hoc  procedure  forms  the  basis  of  the 
present  work. 


II.  FLOW  PROBLEM 

A.  Flow  Equations  and  Boundary  Conditions. 

The  equations  are  stated  and  solved  in  the  inertial  frame.  The  flow  is 
expressed  as  the  sum  of  a  basic  axisymmetric  spin-up  flow  and  a  perturbation:* 


u  =  U  (r,x,f)  -  KQ  u  (r , 9,x,t)  (2.1a) 
v  =  V  (r,x,t)  -  KQ  v  (r.O.x.t)  (2.1b) 
w  =  W  (r,x,f)  -  KQ  w  (r,G,x,t)  (2.1c) 
p  =  P  (r.x.t)  -  KQ  p  (r, 9,x,t) ,  (2. Id) 


where  t  =  $t.  Here  u,  v,  w  are  velocity  components  in  the  radial,  azimuthal, 

★  ★  ★ 

and  axial  directions,  respectively;  U,  V,  W  are  the  corresponding  velocity 

components  of  the  basic  flow  which  is  a  solution  of  the  axisymmetric  Navier- 

★  * 

Stokes  equations.  The  velocity  components  of  the  perturbed  flow  are  u,  v, 

★  ★ 

and  w.  The  quantity  p  is  pressure,  P  is  the  pressure  of  the  basic  flow,  and 
★  • 
p  is  the  perturbation  pressure.  Velocity  is  non-dimensionalized  by  a<}>  and 

pressure  by  pa24>2,  where  p  is  the  density  of  the  liquid. 

The  variables  of  Eq.  (2.1)  are  substituted  into  the  3-D  Navier-Stokes 
equations  (see,  e.g.,  Eqs.  (3.33)  of  Reference  9)  and  only  zeroth  and  first 
order  terms  in  KQ  are  retained.  At  this  point  we  make  several  simplifications 
to  render  the  problem  tractable;  these  steps  are  described  in  detail  in 
Sections  II  and  III  of  Reference  7.  First  the  basic  flow  is  approximated  by 


*  The  negative  siqns  in  Eq.  (2)  were  employed  to  comply  with  the  nomenclature 
of  Reference  3. 

9.  H.  Schlichting ,  Boundary  Layer  Theory ,  McGraw-Hill  Book  Co.,  New  York , 
NY,  1960. 
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the  Wedemeyer  model  for  the  spin-up  flow  (introduced  in  Reference  10  and 
discussed  further  in  Reference  11).  In  this  model  the  flow  field  is  divided 
into  two  regions:  (1)  the  boundary  layers  adjacent  to  the  endwalls,  called  the 
"Ekman  layers,"  and  (2)  the  remainder  of  the  flow  field  called  the  "core,"  in 
which  the  flow  variables  are  designated  U,  V,  W,  P.  (The  need  for  an  addi¬ 
tional  boundary  layer  at  the  sidewall  is  discussed  in  Reference  11.)  These 
satisfy  the  following  set  of  equations  (where  subscripts  denote  partial 


differentiation): 

V£  +  U  (Vp  +  V/r)  =  Re"1  [Vrp  +  (V/r)p]  (2.2a) 

Ux  =  Vx  =  Px  =  °  ("Columnar"  Flow)  (2.2b) 

U  =  fn  (r,V)  (2.2c) 

W  =  -(x/r)  ( rU ) r  (Continuity  Equation)  (2. 2d) 

Pp  =  V2/r  (2.2e) 


Eq.  (2.2b)  indicates  that  V  =  V  (r,  t);  Eq.  (2.2c)  is  the  "Ekman  compatibility 
condition"  given  by  Eqs.  (2.8)  and  (2.9)  of  Reference  7: 


U  =  <  (a/c)  Re"1^  (V-r) 


for  laminar  Ekman  layer 


U  =  -0.035  (a/c)  Re"1^  (r-V)®^  for  turbulent  Ekman  layer; 


<  =  0.5  will  be  used  here.  Boundary  conditions  for  an  impulsive  start  are 

V  (r  -  0)  -  0,  V  (r  -  1)  «  1,  V  (t  =  0)  =  0.  (2.3) 


The  following  linearized  equations  are  obtained  to  describe  viscous 
perturbations  in  the  core  flow  but  not  the  Ekman  layers,  according  to  the 
discussion  preceding  Eqs.  (3.4)  in  Reference  7: 


10.  E.  H.  Wedemeyer ,  "The  Unsteady  Flow  Within  a  Spinning  Cylinder ,"  BRL 
Report  No,  1225,  October  1963,  (AD  431846)  (See  also  Journai_of  Fluid 
Mechanics,  Vol.  20,  Part  3,  1964,  pp,  383-399,) 

11,  R,  Sedney  and  N.  Gerber,  " Viscous  Effects  in  the  Wedemeyer  Model  of 
Spin-Up  From  Rest,"  BRL  Technical  Report  ARBRL-TR-02493 ,  June  1983. 
(AD  A129S06). 
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(2.4a) 


+  (V/r)  uQ  -  2  V  v/r  =  -pp  +  Re"*  (v2u  -  u/r2  -  2  vQ/r2) 

v£  +  (V  +  V/r)  u  +  (V/r)  vQ  =  -pQ/r  +  Re'1  (V2v  -  v/r2  +  2uQ/r2)  (2.4b) 

+  (V/r)  Wg  =  -px  +  Re"1  V2w  (2.4c) 

(ru)r  +  v0  +  r  wx  =  0,  (2.4d) 

where 

V2  H  32/Dr2  +  (1/r)  3/3r  +  (1/r2)  32/302  +  32/3x2. 

The  reasons  that  U  and  W  do  not  appear  in  Eq.  (2.4)  are  discussed  in  Refer¬ 
ence  7. 

The  boundary  conditions  are:  no  flow  through  the  bounding  walls  and  no 
slip  along  them;  i.e., 

u  (r  =  1)  =  w  (r  =  1)  =  0,  v  (r  =  1)  =  1  (side)  (2.5a) 

u  (x  =  ±A)  =  w  (x  =  ±A)  =0,  v  (x  =  ±A)  =  r  (end)  (2.5b) 

where  u,  v,  w  are  non-dimensional  cylindrical  velocity  components  in  the 
r,  6,  x  system;  A  is  the  half-height,  c,  of  the  cylinder  divided  by  a. 

The  boundary  conditions  must  be  transformed  to  the  variables  used  in  Eq. 
(2.4);  a  discussion  of  the  transformations  is  given  in  Appendix  A.  The 
resulting  non-homogeneous  sidewall  conditions  are 

u  (r  =  1)  =  x  Real  [-i  (1-t)  exp  {i  (t^t  -  0)}]  +  0  (KQ)  (2.6a) 

v  (r  =  1)  =  -x  Real  [  (Vp  (r=l)  -t}  exp  (i  (r^t  -  0)}]  +  0  (KQ)  (2.6b) 

w  (r  =  1)  =  Real  [i  (1-t)  exp  (i  (r^t  -  0)}]  +  0  (KQ)  .  (2.6c) 

The  complex  form  is  introduced  here  for  later  convenience. 

★  *  ★ 

The  endwall  boundary  conditions  for  u,  v,  w  cannot  be  stated  precisely. 
The  U,  V,  and  W  are  valid  only  in  the  core,  and  the  Ekman  layer  flow  which 
should  be  included  in  the  unperturbed  solution  is  not  available.  This 
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difficulty  is  responsible  for  the  introduction  of  the  ad  hoc  endwall  boundary 
condition  stated  in  Section  I.  The  motion  of  the  cylinder  is  expected 
nevertheless,  to  determine  the  t,  o  dependence  of  boundary  conditions  so  that 

u  (x  =  ±A)  =  Real  {fn  (±A,r)  exp  [i  (rt  -  e)]},  (2.7) 


★  ★ 

with  similar  expressions  for  v  (x  =  ±A)  and  w  (x  =  ±A).  Axial  boundary 
conditions  will  be  treated  later. 

One  further  approximation  is  needed  to  permit  application  of  a  modal 
analysis,  namely,  the  "quasi -steady"  condition  discussed  in  Section  III  of 
Reference  7.  Generally,  V  does  not  change  appreciably  over  the  time  scale  of 
the  perturbations;  thus  E  can  be  regarded  as  a  parameter  in  the  solution  of 
Eq.  (2.4).  A  possible  exception  to  this  might  exist  for  E  0  and  an 
impulsive  start. 

B.  Modal  Analysis:  Separated  Variable  Solutions 

1.  Equations  and  Boundary  Conditions  for  r,  x  Variation.  It  is  assumed 
that  the  perturbation  can  be  expressed  as  a  superposition  of  modes,  or  a 
triple  Fourier  expansion  in  o,  x,  and  t  with  coefficients  functions  of  r.  It 
is  convenient  to  use  complex  notation  and  express  the  perturbation  as 


u  =  Real  (S  )  *  Real  [£  (r,  x;  E)  exp  (i  (tE  -  e) ) ]  (2.8a) 

v  =  Real  (vc)  =  Real  [v  (r,  x;  E)  exp  {i  (tE  -  9)}]  (2.8b) 

w  =  Real  (5  )  =  Real  [w  (r,  x;  E)  exp  (i  (tE  -  e)}]  (2.8c) 

p  =  Real  (pc)  =  Real  [j)  (r,  x;  t)  exp  {i  (it  -  e)}]  (2.8d) 

*  *  *  * 


where  _u,  v,  w,  and  £  are  complex  functions.  The  functions  uc,  vc,  w£,  and  pc 

are  clearly  also  solutions  of  Eq.  (2.4).  The  t,  0  portion  of  the  boundary 
conditions  in  Eqs.  (2.6)  and  (2.7)  are  satisfied  by  Eq.  (2.8).  Substituting 

it  ★  ★  ★ 

uc»  vc,  wc,  pc  from  Eq.  (2.8)  into  Eq.  (2.4)  yields 


ru^  +  u-  iv  +  rw  =  0 
-r  -  -x 


(2.9a) 


i  (t-V/r)£  -  (2V/rK  =  +  (1/Re)  [urp  +  u^/r  -  2^/r2  +  £xx  +  2i_v/r23  (2.9b) 
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1(T-V/r)v_  +  (Vp+  V/r)  £  =  -ij )/r  +  (1/Re)  [v^r+  Vp/r  - 


(2.9c) 


2v/ r2  +  vxx  -  2iu/r2] 


i(t-V/r)w  =  -£x  +  (1/Re)  [w^  wp/r  -  w/r2  +  wxx].  (2.9d) 

Boundary  conditions  at  the  axis  are 

£  (r  =  0)  -  iv.  (r  *  0)  =  w  (r  =  0)  =  £  (r  =  0)  =  0.  (2.10) 

These  are  kinematic  conditions,  stated  in  Eq.  (4.3)  of  Reference  7,  for  m  =  1. 
Boundary  conditions  at  the  sidewall  are 

£  (r  -  1)  «  -i(l  -  t)x  (2.11a) 

v  (r  =  1)  =  -{Vr(r  =  1)  -t}x  (2.11b) 

w  (r  =  1)  -  i (1  -  t).  (2.11c) 

At  this  point  we  invoke  the  ad  hoc  endwall  boundary  condition,  namely 

w  +  ScE  3w/3x  =  i  (1  -  t)  r  at  x  =  ±A,  (2.12) 


with  the  left-hand  side  identical  in  form  to  that  of  Eq.  (29)  of  Reference  2; 
namely,  the  corrected  homogeneous  endwall  boundary  condition.  The  right-hand 
side  is  the  same  as  that  of  the  uncorrected  inhomogeneous  boundary  condition 
of  Eq.  (A. 10).  The  6cE  is  found  as  follows: 


=  2‘1/2  ReE1/2  (1  -  i)  (3  -  T)1/2 


(2.13a) 


%  = 


2~l/2 


(1  +  i)  (1  +  t)1/2 


(2.13b) 


(2.13c) 
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with  the  restriction  (3  -  t)  >  0.  The  "effective”  Reynolds  number  used  here 
is  given  by 


Re£(E>  = 


V  (r;t)  dr 


a2$/v. 


(2.14) 


For  solid  body  rotation,  with  V  =  r,  the  bracketed  term  is  equal  to  unity  and 
ReE  =  Re.  While  the  fluid  is  spinning  up,  V  <  r,  and  Re£  <  Re.  As  an 

example,  Re^/Re  =  0.76  at  l  =  1000  for  the  case  Re  =  39772,  c/a  -  3.12. 


2.  Form  of  Solution.  As  in  Section  III-B  of  Reference  2,  we  stipulate 
that  the  solution  be  a  TTnear  combination  of  separated  variable  solutions, 
i  .e. , 


u  =  2  Ru(r)  Xu(x) 


v  =  E  Ry(r)  Xy(x) 


w  =  z  Rw(r)  Xw(x) 


£  =  I  R  (r)  X  (x). 


(2.15) 


Substituting  these  into  Eq.  (2.9)  yields  ordinary  differential  equations  for 
the  R's  and  X's.  The  X's  must  satisfy  the  harmonic  equation 


d2X/dx2  +  x2  X  «  0, 


(2.16) 


where  a  can  be  zero  or  finite.  When  x  is  finite,  it  is  determined  from  an 
axial  eigenvalue  problem  and  labeled  with  the  index  k. 

In  this  work  we  choose  to  express  the  solution  as 


KF 


NJ 


£=  u  (r)  x  +  z  u.(r)  sin  A.x  +  Z  d  .  u.(r)  sin  y.) 
o  k  =  1  K  K  j  =  i  J  J  J 


(2.17a) 


A  A 

v  =  v  (r)  x  +  Z  vk(r)  sin  Akx  +  z  d,  v.  (r)  sin  y.x 
-°  k  =  lK  =  J 


(2.17b) 


w  =  w  (r)  -  Z  w.(r)  cos  A.  x  Z  d.  w.(r)  cos  u.x 
-  o  k  =  i  k  k  j  =  i  J  J  J 


(2.17c) 


£  =  P0(r)  x  +  E  pk(r)  sin  A^x  +  E  d^  pj(r)  sin  y^x, 
k  -  1  j  *  1 


(2.1 7d ) 
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where  the  x^'s  form  a  denumerable  set  of  solutions  to  the  functional  equation 
(See  Eq.  (31)  of  Reference  2) 

cos  X^A  Xk  Sc£  sin  =  0.  (2.18) 

For  | 6ce | /A<<1 , 

\  =  (ku)/[2  (A  -  6c E ) ] .  k  odd  (2.19) 

A 

Thus,  at  the  endwall  each  wk(r)  cos  xkx  term  in  Eq.  (2.17c)  satisfies  Eq. 
(2.12)  with  the  right-hand  side  replaced  by  zero. 

The  y-'s  are  eigenvalues  of  the  radial  differential  equations  resulting 

J 

from  separation  of  variables,  with  homogeneous  sidewall  boundary  conditions: 

u.  (r  =  1)  =  v.  (r  =  1)  =  w.(r  =  1)  =  0.  (2.20) 

J  J  J 

The  parameter  t  enters  here  in  the  same  way  as  the  unknown  eigenvalue  C  in 
Reference  7.  In  the  present  situation  t  is  known  and  y.  is  the  eigenvalue  to 

J 

be  determined.  The  dj's  are  constants  to  be  determined  so  as  to  minimize  the 

error  in  satisfying  the  endwall  boundary  condition,  Eq.  (2.12).  The  homogene¬ 
ous  boundary  conditions,  Eqs.  (2.10)  and  (2.20),  do  not  produce  unique  eigen¬ 
functions  u j ,  Vj,  w j ,  and  p^;  these  are  known  to  within  a  common  constant. 

The  functional  values  obtained  from  the  integration  are  a  result  of  the  values 

assigned  the  non-homogeneous  quantities  at  r  =  0.  The  subsequent  determina¬ 
tion  of  the  dj's,  however,  renders  the  flow  solution  unique. 

A  A  A  A 

The  uQ,  vQ,  wQ,  pQ,  corresponding  to  X  =  0,  satisfy  the  sidewall  condi¬ 
tions  (the  same  as  for  the  particular  solution  in  Reference  2) 

uQ  (r  =  1)  =  -i  (1-t)2/(1+t)  (2.21a) 

vQ  (r  =  1)  =  — C(l-r)2/ (1+t)]  (2.21b) 

wq  (r  =  1)  =  i  (1-t).  (2.21c) 

A 

This  set  of  conditions  was  chosen  so  that  wQ  satisfied  the  total  boundary 
condition  at  r  =  1,  Eq.  (2.11c);  better  conditions  on  uQ(l)  and  vQ(l)  might  be 
selected,  and  alternate  choices  are  being  considered.  The  sidewall  boundary 
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conditions  for  second  terms  of  the  solution,  Eq.  (2.17),  are  obtained  by 

A  A 

subtracting  the  expressions  of  Eq.  (2.21)  (times  x  for  uQ  and  vQ)  from  those 
of  Eq.  (2.11) : 

Z  uk  (r  =  1)  sin  Akx  =  -i  [2t  (1-t)/(1+t)]  x  (2.22a) 

£  vk  (r  =  1)  sin  ^x  =  [ {(1-t)2/(1+t) }  -  Vp  (r  -  1)  +  t]  x  (2.22b) 

Z  wk  (r  =  1)  cos  Akx  =  0.  (2.22c) 

For  solid  body  rotation  the  dj's  vanish  and  the  remainder  of  the  solution,  Eq. 
(2.17),  is  the  same  as  that  in  Reference  2. 

3.  Radial  Variation  Problem.  When  each  individual  term  of  Eq.  (2.17) 
is  substituted  into  tq.  (2.^),  a  system  of  linear  ordinary  equations  in  r  of 
the  fol lowing  form,  omitting  subscripts  k  on  the  dependent  variables,  is 
obtained,  (where  '  =  d/dr) : 

AAA  A 

r  u'  +  u  -  iv  +  Ak  rw  =  0  (2.23a) 

Re"*  u"  +  (Re  r)'1  u‘  +  [i  (V/r  -  t)  -  Re'1  (2/r2  +  Ak2)]  u  + 

(2.23b) 

[2  V/r  +  2i  (Re  r2)"1]  v  =  p' 

Re'1  v"  +  (Re  r)"1  v'  +  [i  (V/r  -  T)  -Re'1  (2/r2  +  Ak2)]  v  - 

(2.23c) 

[V/r  +  Vp  +  2i  (Re  r2)'1]  u  =  -i  p/r 

Re'x  w"  +  (Re  r)'1  w'  +  [i  (V/r  -  t)  -  Re"1  (1/r2  +  Xk2)]  w  = 

(2.23d) 

'*k  P  +  ep  P* 
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The  quantity  Rp  is  equal  to  1  when  x^  =  0;  otherwise,  it  is  zero.  The  same 
equations  apply  to  u^,  Vj,  w^,  p^,  with  Xk  replaced  by  uj. 

These  equations  are  converted  to  canonical  form  in  order  to  be  integrated 
numerically;  i.e., 

y,  =  dy^/dr  -  (r,  yp  y2»»»«yg)»  '  =  1,2,. ..,6, 


where 


A  A 

y^  =  u  (or  u)  y4  =  w  (or  w) 

A  A  A  | 

y2  =  u  -  iv  (or  u  -  iv)  y&  =  w'  (or  w' )  \  (2.24) 

y3  =  v'  (or  v')  y6  *  P  (or  p).  ' 

After  the  required  manipulations  are  performed,  the  following  sixth  order 

A  A  A  A 

system  is  obtained  for  the  u,  v,  w,  p  case: 


V  =  (y2/r>  -  \  y4 

y2‘  =  -  (y2/r)  -1'  y3  -  xk  y4 

y3‘  »  [Re  (V/r  +  Vp)  +  2i/r2]  yx  +  i  (B  +  1/r2)  (y2  -  yx)  - 

y3/r  -  i  (Re/r)  yg 

y4*  =  y5 

V  =  B  *4  ‘  V  '  Xk  Re  >6  +  ep  Rc  y6 

y6‘  =-Byj/Re  ♦  ( i /Re)  [2  Re  V/r  +  i/r2]  (y2  -  yj)  ♦ 

i  y3/(r  Re)  -  xk  y5/Re, 


(2.25a) 

(2.25b) 

(2.25c) 

(2.25d) 

(2.25e) 

(2.25f) 
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where 


B  =  (1/r2)  +  X,^2  +  i  Re  [t  -  V/r], 


(2.26) 


The  X^'s  are  replaced  by  p^'s  for  the  u,  v,  w,  p  case. 


4.  Boundary  Conditions  for  Radial  Equations.  There  are  three  boundary 
conditions  at  r  =  0  and  three  at  r  =  1.  As  a  consequence  of  Eqs.  (2.10), 
(2.17),  and  (2.24),  conditions  at  r  =  0  are 


y2k(o)  =  y4k(o)  •  y6k(o)  =  o 

y2J(°)  -  y4j<°)  ■  y6j«»  -  o. 


In  Eq.  (2.22)  the  function  x  can  be  expanded  in  the  series 


x  *  I  sin  x^x, 
k 


(2.28) 


where  the  coefficients,  given  by  Eq.  (40)  in  Reference  2,  are 


(2/  X^2)  [1  +  (Xk  6ce)2]  sin  xk  A 
k  A  [1  +  (Xk  6Ce)*]  -  6CE 

(2.29) 

These  are  the  coefficients  in  a  series  of  biorthogonal  functions  determined  by 
solving  a  non-self-adjoint  system,  Eq.  (2.16).  By  applying  Eqs.  (2.21), 
(2.22),  (2.24),  and  (2.29),  we  obtain  the  following  sidewall  conditions  for 

the  uk,  vk,  wk,  pk  case: 

ylk(l)  =  -i  (l-cp)  [2t  (1-t)/(1+t)]  bk  -i  £p  [(1-t)2/(1+t)] 

(2.30a) 

y2k(l)  =  -i  (1-Ep)  [1  -  Vr(l)]  bk 

(2.30b) 

y4k(!)  s  i  ep  U-t). 

(2.30c) 

For  the  eigenvalue  problem  solutions,  u.,  v.,  w.,  p.  case, 

J  J  J  J 
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(2.30d) 


ylj(1)  "  y2j^  =  y4j^  =  °* 


5.  Endwall  Boundary  Condition.  When  we  substitute  the  w  solution,  Eq. 
(2.17c)  evaluated  at  x  *  ±  A,  into  Eq.  (2.12)  and  apply  Eq.  (2.18)  to  the 


cos 


x  terms,  we  arrive  at  the  following  endwall  boundary  condition: 


NJ 


F  (r,  ±A)  =  ^  (r)  -  l  d.  w.(r)  (cos  u-  A  +  6cF  u,  sin  p.  A)  =  0, 

j  =  J  J  J  J  J 


where 


and 


*  (r)  *  w Q(r)  -  i  (1-t)  r  =  ^  +  i  ^ 


F  =  ER  +  l  Fr 


(2.31) 


(2.32) 

(2.33) 


As  discussed  in  Section  I,  Eq.  (2.31)  is  not  satisfied  exactly;  instead  we 
determine  a  set  of  dj's  which  will  minimize 

9(dj»  d2»*..d^j)  =  ^  (Fr  +Fj)dr.  (2.34) 

A  derivation  of  the  formulas  for  the  dj's  is  given  in  Appendix  B.  Thus,  the 

ad  hoc  endwall  boundary  condition  is  satisfied  approximately,  to  within  the 
accuracy  of  a  least  squares  fit.  An  estimate  of  the  relative  error  is  given 
by  the  quantity 

Er  =  [g/  (  UR2  +  *j2)  dr]1/2.  (2.35) 

Error  estimates  have  been  computed,  but  we  defer  the  presentation  of  results 
to  a  later  report. 

C.  Operational  Procedures. 

The  procedure  for  solving  the  radial  ordinary  differential  equations,  Eq. 

(2.25),  is  discussed  in  detail  in  Reference  1  for  solid-body  rotation  and  in 
Reference  7  for  spin-up.  Numerical  integration  must  begin  at  a  small 

finite  value,  r  =  e0,  where  three  independent  power-series  solutions  are 

evaluated.  Due  to  the  stiffness  of  the  equations  at  large  Re,  the  Runge-Kutta 
integration  must  be  combined  with  orthonormalization  of  the  three  solutions  in 
order  to  maintain  their  linear  independence  and  prevent  runaway  amplitude 
growth.  Both  homogeneous  and  inhomogeneous  sidewall  conditions  are 
encountered  in  this  study.  Reference  1  provides  details  of  the  application  of 
the  homogeneous  conditions,  and  Reference  2  describes  the  treatment  of  the 
non-homogeneous  conditions.  The  iterative  procedure  for  computing  eigenvalues 
is  described  in  Section  V  of  Reference  7. 
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Because  the  calculation  of  the  pj  eigenvalues  forms  part  of  our  opera¬ 
tions,  it  is  most  feasible  to  begin  the  determination  of  a  spin-up  moment 
history  at  a  very  large  time,  near  solid  body  rotation,  where  initial  guesses 
for  the  eigenvalues  are  readily  available.  These  are  found  in  tables  (see, 
e.g.,  p.  31  of  Reference  3).  The  operation  procedes  step-by-step  from  larger 
to  smaller  values  of  time,  with  initial  guesses  for  the  eigenvalues  at  each 
time  step  obtained  from  those  at  previous  times. 


For  large  t  the  p.'s  can  be  ordered  (j  =  1,2,  etc)  either  by  the  real 

J 


parts  of  the  p.'s  or  the  number  of  zeros  in  their  corresponding  eigenfunc- 

J 

tions.  The  mode  is  then  identified  by  the  value  of  j.  However,  as  t  becomes 
small,  some  of  the  Real  (p.)  vs  t  curves  often  intersect,  rendering  the  ordering 

J 

by  Real  (p.)  inapplicable.  At  these  times  identification  of  the  modes  can 
J 


become  ambiguous.  What  frequently  happens  is  that  the  iterative  process 
converges  to  the  solution  of  a  mode  other  than  the  one  sought  because  the 
initial  estimate  was  not  close  enough  to  the  desired  eigenvalue.  Opera¬ 
tionally,  this  can  lead  to  a  particular  mode  appearing  more  than  once  in  the 

expansion  Ed.  u.  (r)  sin  p.  x,  etc  (Eq.  2.17a)).  We  can  remove  one  of  the 
3  3  3 


duplicate  terms  from  the  series,  but  we  are  then  left  with  fewer  functions  for 
approximating  the  endwall  boundary  condition.  As  the  calculation  proceeds  for 
decreasing  time,  the  number  of  eigenfunctions  decreases  by  one.  The  error, 
Er,  then  increases. 


Another  possible  source  of  difficulty  at  early  time  occurs  when 
M  =  t-  V/r  =  0,  i.e.,  when  the  nutational  frequency  equals  the  frequency  of 
the  circumferential  motion  of  the  fluid,  indicating  a  resonance;  see  Refer¬ 
ences  7  and  12.  The  value  of  r  for  which  M  vanishes  is  called  the  "critical 
level"  and  is  denoted  by  rc;  the  neighboring  region  is  called  the  "critical 

layer."  The  critical  level  begins  at  r  =  1  for  t  =  0  and  moves  inward  with  in¬ 
creasing  time  until  it  reaches  the  axis  at  a  finite  time  and  then  disappears. 
When  inviscid  perturbations  are  treated  for  constant  yaw,  as  in  Reference  8, 
the  flow  equations  become  singular  at  r  =  rc.  The  critical  layer  is  discussed 

in  References  7  and  12  for  viscous  perturbations.  Figure  4c  in  Reference  7 
illustrates  in  this  situation  that  high-frequency  large  amplitude  oscillations 
in  the  flow  solution  occur  about  r  =  rc  under  certain  conditions;  the  numeri¬ 
cal  process  described  in  Reference  7  is  able  to  compute  these  oscillations. 
The  number  of  zeros  becomes  large  and  identification  of  a  mode  is  again  made 
ambi guous. 


III.  LIQUID  PRESSURE  MOMENT 
A.  Pressure  and  Moment  Formulas. 

According  to  Eqs.  (2. Id)  and  (2.8d)  the  pressure  is 


12.  R.  Sedney  and  N.  Gerber,  "Numerical  Study  of  the  Critical  Layer  in  a 
Rotating  Fluid,"  AIAA  22nd  Aerospace  Sciences  Meeting,  Reno,  Nevada, 
January  1984,  AIAA  Paper  No.  84-0342. 
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(3.1) 


p  =  P  -  Kq  [-gj  sin  ( x*t-o)  +  gR  cos  (x*t-o)]  +  0  (k/), 

where 


P  (r»x)  0R  i  pj.  (3.2) 

Omitting  the  despin  torque,  one  can  express  the  moment  produced  by  the  liquid 
on  the  spinning  and  nutating  shell  as  the  complex  quantity  M^~jp  +  i  M^L^p. 

A  complex  moment  coefficient,  C/)M\P,  is  introduced  as  in  Eqs.  (5)  and  (6)  of 
Reference  2:  {  ’ 

M(LY)P  +  1  M(LZ)P  =  (2"0a2c)  a^2  T  C(LM)P  K1  e  ^ 

(3.3) 

C(LM)P  =  C(LSM)P  +  1  C(LIM)P* 

The  i_sm) P  an<:1  C(|_IM)P  represent  the  moments  acting  to  change  the  yaw  angle 
and  the  nutation  rate,  respectively. 

The  moment  will  be  evaluated  about  the  center  of  gravity  of  the  projec¬ 
tile  in  the  x,  y,  z  system.  Details  need  be  shown  for  only  one  component, 
say  M(l£)P’  because  °f  axisymmetry  in  the  transverse  motion.  Let 

M(LZ)P  =  M(LZS)P  +  M(LZT)P  +  M(LZB)P*  (3*4) 

where  the  three  terms  on  the  right-hand  side  denote  the  moments  on  the  side, 
top,  and  bottom  walls,  respectively;  as  given  in  Eq.  (53)  of  Reference  2  they 
are 


'fA 

r 2" . ■ 

M(LZS)P 

=  (pa5*2) 

\  p  (r  =  1)  x  cos  o  de  dx 

J  -A 

° 

[V  i 

.  2u 

M(LZT)P 

=  -(pa5*2) 

1 

l  p  (x  =  A)  r2  cos  0  de  dr 

.Jo  • 

o 

■f1  f 

2  it 

M(LZB)P 

=  (pa5*2) 

p  (x  =  -A)  r  cos  o  de  dr 

J  0  J 

0 

(3.5a) 


(3.5b) 


(3.5c) 
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The  non-dimensional  cylinder  radius  occurring  in  the  integrand  of  Eq.  (3.5a) 
is  equal  to  1. 

B.  Sidewall  Moment. 

The  pressure  must  be  evaluated  at  r  =  1.  The  bracketed  term  in  Eq.  (3.1) 
may  be  evaluated  at  r  =  1  without  changing  the  order  of  the  approximation. 
According  to  Eq.  (A4) 

P  (r  *  1)  *  P  (r  *  1)  -  Kqx  cos  (x$t-e)  (aP/3r)r  _  j. 

P  (r  =  1)  is  a  constant,  designated  by  Pj.  From  Eq.  (2.2e),  (3P/3r)r  _  ^  = 
(V2/r)r  _  j  =  1.  Hence 

P  (r  =  1)  =  P  (r  =  1)  -  KQ  x  cos  (x$t-e).  (3.6) 

Then,  from  Eqs.  (2. Id),  (2.8d),  and  (3.6),  with  P  replacing 


p  (r  «  1)  »  Pj  -  KQ  {[x  +  pR  (l,x)]  cos  (x$t-e)  - 
Pj  (l,x)  sin  ( x^t-O) }  +  0  (Kq2). 


(3.7) 


When  Eq.  (3.7)  is  substituted  into  Eq.  (3.5a),  the  term  makes  no 

contribution  to  the  integral,  leaving  only  KQ  terms;  thus  r,  e,  x  can  be 

replaced  by  r,  o,  x  without  changing  the  order  of  approximation.  Then  Eq. 
(3.5a)  reduces  to 


M(L2s)p/(K0  pa5*2)  =  -[2*  Real  (Is)]  sin  x$t  - 

[2n  Imag  (1^)]  cos  x$t 


where 


C  ^ 

Is  =  i  A3/3  +  i  l  x  e  (1.x)  dx. 

'  0 

Upon  substitution  of  Eq.  (2.1 7d ) .  Eq.  (3.9)  becomes 


(3.8) 


(3.9) 
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(3.10) 


I  =  1  A3  Cp  ( 1)  +  1]  +  i  E  p  (1)[_L  sin  x.  A  -  —  A  cos  X.A]  + 
b  3  0  k  =  1  k  x2  k  XR  k 


i  E  d.  p.(l)  [-J—  sin  p.A  -  —  A  cos  u.A]. 

j  =  1  J  J  g.2  J  u.  J 


C.  Endwall  Moments. 

The  term  P(r)  of  Eq.  (3.1)  must  be  evaluated  as  a  function  of  r,  e  at  x  = 
A  in  the  moment  integral  of  Eq.  (3.5b)  for  the  top  endwall.  As  demonstrated 
in  Appendix  A,  Eq.  (A7), 


P(r)  =  P(r)  -  K0  A  (V2/r)  cos  (T*t-o)  +  0  (KQZ) .  (3.11) 


The  P(r)  term  makes  no  contribution  to  the  integral,  leaving  only  terms  of 
0(Ko);  then  the  r,  e,  x  may  be  replaced  by  r,  6,  x  without  changing  the  order 


of  approximation.  The  remaining  part  of  the  pressure  is  an  odd  function  of  x; 
thus,  the  total  endwall  moment,  m(qe)p  =  m(lZT)P  *  M(LZB)P  =  2  M(LZT)P* 


Eqs.  (3.1)  and  (3.11)  are  substituted  into  Eq.  (3.5b)  to  yield 
M( lze ) p/ ( p3 5<^2)  ~  2l*  ko  (t Imag  (1^)3  cos  t<j>t  + 

[Real  ( 1^) ]  sin  t$t} 

where  ^ 

Ie  =  i  j  r2  [p  (r,  x  =  A)  +  A  V2/r]  dr. 

Jo 

Substituting  Eq.  (2.1 7d )  into  Eq.  (3.13)  gives 

f1 

IE  =  I  r2  [A  pQ  (r)  +  E  pk  (r)  sin  Xk  A  + 

Jo  k  =  1 

NJ 

E  d.  p.  (r)  sin  p.  A  +  i  A  V2/r]  dr. 
j  =  1  J  J  J 


(3.12) 


(3.13) 


(3.14) 
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D.  Moment  Coefficients, 

The  moments  of  Eqs.  (3.8)  and  (3.12)  have  the  form 


M(LZS)P  =  ”Ko  (M1S  sin  +  M2S  cos  *♦*) 

(3.15) 

M(L2£)p  =  -K0  (M^  sin  x$t  +  M££  cos  t <j>t )  (pa^<}>^), 

where 


'lS  =  ^ea^  (Is)- 


Mgj  =  2n  Imag  (I<j) 


=  “2k  Real  (I^),  ^2E  =  Imag  (I^)< 


The  total  moment  is 


'(LZ)P  =  ^(LZS)P  *  ^(LZE)P  =  -Ko  (M^sin  x$t  +  M2  cos  x4>t),  (3.16) 


where 


Mi  =  M1S  +  M1E, 


M2  =  M2S  +  M2E* 


(3.17) 


Our  computational  results  will  be  exhibited  in  terms  of  the  moment 
coefficient,  C(Lm)p*  defined  in  Eq.  (3.3).  By  Eqs.  (3.3),  (3.16),  and  (3.17), 

C(lsm)P  =  ~M^/C2kx  c/a],  ^(LIM)P  =  c/a].  (3.18) 

Our  primary  interest  lies  in  C(i_$M)p,  which  represents  overturning  moment,  and 
we  shall  present  time  histories  of  this  quantity. 


IV.  RESULTS 

This  report  does  not  include  the  contribution  to  the  moment  from  the 
shear  forces  which,  as  shown  in  References  3  and  4,  can  be  significant  at 
lower  Re;  shear  contribution  will  be  treated  in  a  subsequent  work.  For  the 
present  we  concentrate  on  the  analysis  of  the  flow  problem  and  procedures  for 
solving  it,  which  are  needed  for  both  pressure  and  shear  contributions.  The 
results  given  here  illustrate  the  procedure  and  are  compared  with  those  of 
inviscid  perturbation  calculations  by  the  method  of  Murphy.8  Side  pressure 
moment  coefficient  histories  are  presented  for  four  combinations  of  Re  and 
c/a. 
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histories  are  shown  for  five  nutational 


In  Case  A,  Figure  2,  C(lsm)p 

frequencies:  x  =  0.04,  0.05,  0.09,  0.12,  and  0.14.  The  curves  have  very 
pronounced  peaks  for  t  =  0.05,  0.09,  0.12,  and  0.14.  These  peaks  (except  for 
the  right-hand  peaks  on  the  t  =  0.12  and  t  =  0.14  curves)  occur  at  times  close 
to  those  at  which  the  perturbing  motion  is  in  resonance  with  the  wave  motion 
of  the  n  =  1,  k  =  3  natural  oscillation  mode  of  the  liquid  in  the  cylinder.* 
The  right-hand  peaks  on  the  r  =  0.12  and  x  =  0.14  curves  indicate  resonance 
with  the  n  =  2,  k  =  5  mode  of  the  spinning  fluid.  The  eigenf requency  of  the  n 
=  1,  k  =  3  mode  is  approximately  Cp  =  0.04  for  solid  body  rotation; 

consequently  the  moment  will  not  damp  out  with  increasing  time  for  x  =  0.04. 
Outputs  of  the  inviscid  perturbation  calculations  (method  of  Reference  8)  are 
also  presented  for  x  =  0.04,  0.05,  and  0.09.  The  predictions  of  times  of  peak 
moments  differ  by  less  than  0.07  s.  between  the  two  methods.  For  x  =  .12  and 
.14  only  the  present  results  are  shown.  Reference  8  cannot  compute  these  cases 
because  of  the  presence  of  the  critical  layer. 

The  results  for  Case  B  are  shown  in  Figure  3.  The  spin-up  time  (defined 
in  Section  I  of  Reference  11)  is  smaller  here  than  in  Case  A;  therefore,  the 

interval  of  interest  occurs  for  smaller  £.  The  results  are  qualitatively  the 

same,  however.  The  maximum  of  C(lsm)p  occurs  near  the  time  when  Cp  =  x  =  0.15 

for  the  k  =  3,  n  =  1  mode.  The  results  from  the  method  of  Reference  8  are 

also  shown  in  Figure  3.  These  disagree  with  the  present  results  at  t  =  200. 

For  smaller  £  that  method  breaks  down  because  of  the  existence  of  the  critical 
layer.  Our  calculation  was  not  carried  out  for  smaller  times  because  of 
difficulties,  described  in  Section  II. C,  which  were  encountered  in  evaluating 
the  ,j.'s. 

J 

In  Case  C  the  high  Reynolds  number  requires  a  turbulent  Ekman  layer.  The 

aspect-ratio  and  nutational  frequency  are  chosen  so  that  Cp(t)  =  x  at  £  = 

38,500  for  the  k  =  5,  n  =  1  mode.  Both  methods  exhibit  the  occurrence  of 
resonance  near  this  time.  The  critical  layer  occurs  beyond  the  range  of  these 
curves  in  this  graph. 

The  only  difference  between  the  parameters  of  Case  C  and  Case  D,  Figure 
5,  is  a  2.55%  increase  in  aspect  ratio,  yet  the  moment  coefficient  histories 
are  radically  different.  Sensitivity  to  aspect  ratio  has  been  noted  previ¬ 
ously  in  yaw  growth  rates  of  projectiles  containing  liquid  payloads  in  solid- 
body  rotation;  see,  e.g.,  Figures  10  and  12  in  Reference  4.  Both  methods  show 
yaw  damping  moments  at  this  stage  of  the  spin-up,  with  maximum  damping  near 
the  time  when  Cp  =  x;  however,  the  amplitudes  differ  widely.  Furthermore,  the 

curve  of  Murphy's  method  exhibits  a  period  of  sharply  rising  and  falling 
overturning  moment  which  the  present  method  does  not.  The  results  of  Case  D 
are  not  yet  thoroughly  understood.  The  following  comments  may  be  relevant  to 


AThe  reader  is  referred  to  References  1  and  7  for  description  of  the  oscilla¬ 
tion  modes  of  a  liquid  cylinder.  In  this  study  spin-up  eigenvalues  were  com¬ 
puted  without  endwall  boundary  condition  correction. 
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the  situation.  Because  the  cylinder  is  executing  forced  rather  than  free 
angular  motion,  the  liquid  may  under  certain  conditions  exert  a  yaw  damping 
moment.  The  current  treatment  of  the  endwall  boundary  condition  would  lead 
one  to  infer  that  the  theory  grows  less  valid  with  decreasing  time.  In  the 
inviscid  perturbation  method  the  results  may  deteriorate,  when  operating 
backwards  in  time,  before  the  critical  layer  first  appears. 

For  each  perturbing  frequency  there  is  a  time  t  =  tM  at  which  the  moment 

coefficient  has  a  maximum,  C^LSMjPM.  Conversely,  at  each  time  there  is  a 

perturbing  frequency  (t)  for  which  C(LSM)p  =  C(LSM)PM’  For  Case  A  Figure 

6  exhibits  the  upper  limit  to  the  side  moment  coefficient  that  can  be  experi¬ 
enced  during  spin-up  compared  to  that  which  can  be  experienced  for  fully  spun- 
up  basic  flow.  The  moment  coefficient  ratio  decreases  rapidly  with  decreasing 
time  in  the  present  calculation,  but  the  inviscid  perturbation  result  indi¬ 
cates  a  very  slowly  changing  ratio. 


V.  SUMMARY 

The  interaction  between  the  flow  of  a  liquid  payload  and  the  motion  of 
the  spinning  and  nutating  shell  containing  it  complicates  the  prediction  of 

the  angular  motion  of  the  projectile  and  the  forces  on  it.  This  paper  repre¬ 
sents  a  first  effort  using  linearized  fully  viscous  equations  to  determine 
side  pressure  moments  exerted  by  the  liquid  on  the  container  during  the  liquid 
spin-up  process.  The  heuristic  approach  of  Murphy8  is  used  which  makes 

crucial  assumptions  regarding  time-dependence  of  the  flow  and  the  endwall 

boundary  condition.  The  statement  of  the  problem  specifies  the  angular  motion 
of  the  projectile;  thus,  the  present  treatment  cannot  simulate  actual  flight, 
although  it  can  simulate  performable  gyroscope  experiments. 

Calculations  were  made  for  four  cases.  It  is  seen  that  yaw  damping 
moments  exist  under  certain  circumstances.  For  three  cases  good  qualitative 
agreement  is  obtained  with  the  results  of  inviscid  perturbation  calcula¬ 
tions.  Side-moment  coefficient  histories  are  bell-shaped  curves  with  peaks 

occurring  approximately  at  times  when  the  nutational  frequency  coincides  with 
an  eigenf requency  of  the  liquid.  (For  the  solid-body  rotation  resonance  case 
in  Figure  2  the  curve  has  a  flat  top.)  In  the  fourth  case  the  discrepancy 

between  the  results  of  the  two  methods  is  large,  though  both  methods  predict 

sharp  peaks  of  damping  moment  coefficient. 

The  next  step  in  this  investigation  is  to  compute  shear  force  on  the 

cylinder  walls  and  then  obtain  the  shear  moment  on  the  projectile.  It  is 

expected  to  have  a  noticeable  effect  for  smaller  Reynolds  numbers.  Inclusion 
of  the  Ekman  layers  in  the  basic  flow  is  necessary  in  order  to  have  a  rational 
approximation  to  the  solution. 
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Figure  1.  Diagrams  of  Coordinates  and  Cylinder. 
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Figure  2.  Side  Moment  Coefficient  Histories  for  Five  Nutational  Frequencies, 
Case  A:  Re  =  39772,  c/a  =  3.12,  $  =  754  rad/s. 
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Figure  3.  Side  Moment  Coefficient  History,  Case  B:  Re  =  4974, 
c/a  *  3.30,  *  8937  rad/s,  x  *  0.15. 


Figure  4.  Side  Moment  Coefficient  History,  Case  C:  Re  =  1.85  x  106, 
c/a  *  5.100,  {  =  641  rad/s,  x  =  0.025. 
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C(LSM)P 


Figure  5.  Side  Moment  Coefficient  History,  Case  D:  Re  =  1.85  x  106, 
c/a  =  4.973,  $  =  641  rad/s,  t  =  0.025. 


C(LSM)PM 
C(LSM)PM  ^ 


Figure  6.  Moment  Coefficient  Ratio  vs  $t  for  Case  A: 
Re  =  39772,  c/a  =  3.12. 
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APPENDIX  A:  EVALUATION  OF  VARIABLES  ON  BOUNDARIES 


The  relationship  between  the  inertial  and  nutating  sets  of  cylindrical 
coordinates  is  (Eq.  (9)  in  Reference  2) 


r  =  r  -  x  cos  (x$t-e)  +  0(Kq2) 

9  =  0  -  Kj  (x/r)  sin  (x$t-e)  +  0(KQ2) 

x  =  x  +  Kj  r  cos  (x$t-e)  +  0(Kq2)  . 


(A. la) 
(A. lb) 
(A.lc) 


Recalling  Eq.  (2.1),  velocities  are 


u  =  dr/df  = 


★  * 

U  -  K0u, 


v  =  r  de/df  * 


*  * 
V  -  K  v 
o 


w  =  dx/dt  = 


★  ★ 
W  -  Kqw 


(A. 2) 


u  =  dr/dt. 


v  =  r  de/df. 


w  =  dx/dE. 


The  velocity  transformation  between  the  earth-fixed  and  aerobal 1 istic  systems 

is  obtained  by  differentiating  Eq.  (A.l)  with  respect  to  time  and  substituting 

the  expressions  of  Eq.  (A. 2)  for  the  derivatives.  Here  we  make  an  additional 

approximation  in  order  to  keep  the  expressions  tractable.  According  to  the 

-1/2  * 

Wedemeyer  model ,  U  and  W  are  both  of  order  Re  ,  we  assume  the  same  for  U 

*  -1/2  2  1/2  2 

and  W.  We  shall  neglect  KQ  Re  '  and  KQ  Re  '  terms  in  addition  to  KQ 

terms.  We  then  obtain  for  e  =  0,  i.e.,  =  Kq, 


u=U-Kou+Kqx  (V/r  -  x)  sin  (x^t-0)  (A. 3a) 

v  =  V  -  Kqv  +  Kq  x  t  cos  (T4>t-e)  (A. 3b) 

w  =  W  -  Kqw  -  KQr  (V/r-x)  sin  (x4>t-e).  (A. 3c) 

A  variable  may  be  evaluated  at  r  =  1  by  an  expansion  about  r  =  1,  thus 


h  (r  =  1)  =  h  (r  =  1)  +  (3h/ar)r  =  j  (r“  =  1}  -  1)  +  0(KQ2). 
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By  application  of  Eq.  (A. la) 


h  (r  *  1)  =  h  (r  =  1)  -  KQx  cos  (i^t-0)  (3h/3r)r  _  j  +  0( KQ2) .  (A. 4) 


Equations  (A. 3)  are  now  evaluated  at  r  =  1.  Recalling  that  U,  l5  =  0(Re~^2) 
and  U  (r  =  1)  =  W  (r  =  1)  =  0,  we  apply  Eq.  (A. 4)  to  Eq.  (A. 3)  to  obtain  at 
r  =  1: 


0  =  Kq  [-u  +  x  ( 1  - t)  sin  (t-jit-e)]  +  0  (Kq2,  Re"1/2  KQ)  (A. 5a) 
1  ■  1  +  Kq  [-v  -  x  (Vr(r  =  1)  -  t)  cos  (4t-0)]  (A. 5b) 
0  s  Kq  [-w  -  (1-t)  sin  (x^t-e)].  (A. 5c) 


Finally,  replacing  0  in  Eq.  (A. 5b)  by  V  of  the  Wedemeyer  model  we  obtain  the 
boundary  conditions  of  Eq.  (2.6).* 


To  evaluate  P  on  the  endwall,  one  applies  the  expansion 
P  (x  =  A,  r)  =  P  (x  =  A,  r)  +  (x“  =  A  -  A)  (3P/3x)x  _  A  + 


(r  -  r)  (3P/3r)p  =  p  +  0(Kq2). 


(A.6) 


*Eqs.  (A2)t  ( A4 ) ,  and  (AS)  in  Reference  2  contain  errors.  The  signs  directly 

a  *  * 

following  u,  v ,  and  u  in  the  right-hand  sides  of  Eqs.  ( A2)  should  all  be 
positive.  The  minus  signs  directly  following  the  equals  signs  in  Eqs.  (A4) 

and  (A5)  should  all  be  changed  to  plus  signs.  Also  v  (r  =  1)  in  the  last 

equation  of  Eqs.  (A4)  should  be  replaced  by  w  (r  =  1). 
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In  the  Wedemeyer  model  (3P/3x)  =0;  (r  -  r)  is  given  by  Eq.  (A. la).  Thus 


P  (r,  x  =  A)  =  P  (r)  +  KQ  A  (V2/r)  cos  (x$t-e)  +  0(KQ2).  (A. 7) 

The  condition  of  zero  normal  flow  at  the  endwall  is  obtained  by  evaluat- 

★  ★  - 

ing  Eq.  (A. 3c)  at  x  =  A,  where  W  »  0  and  V/r  *  1.  The  term  w  (x  =  A)  is  zero 
in  our  approximation.  This  is  shown  in  the  expansion 

w  (x  =  A)  =  w  (x  =  A)  +  [x  (x  =  A)  -  A]  (3w/3x)“  _  A  +  0(KQ2).  (A. 8) 

By  Eq.  (2.5b),  w  (x"  =  A)  =  0;  by  Eq.  (A.lc).  [x  (x  =  A)  -  A]  =  0(Kq);  by  Eq. 

(A. 3c),  w  =  W  +  0(K  ).  In  the  Wedemeyer  model  W  and  its  derivatives  are  of 

-1/2  0  *  ~  . 
order  Re  ;  the  same  is  assumed  for  W  and  its  derivatives.  Then  w  (x  =  A)  = 

-1/2 

0(Kq  Re  '  )  and  is  neglected,  leaving 

w  (x  =  A)  =  -(1-t)  r  sin  (x$t-e),  (A. 9) 

and 

w  (x  =  A)  =  i  (1-x)  r.  (A. 10) 
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APPENDIX  B 

FORMULAS  FOR  dj's  IN  EQ.  (2.17) 
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APPENDIX  B:  FORMULAS  FOR  dj's  IN  EQ.  (2.17) 


We  define  ej  by 


ej  5  eRj  +  *  eIj  "  (cos  “j  A  +  {CE  “j  Si"  "j  A)  dj 


(B.l) 


Eq.  (2.31)  is  now  rewritten: 


NJ 


FR  +  1  FI 


*  (r)  -  l  [eRj.  +  i  e^]  [wRj.  (r)  +  i  w^  (r)],  (B.2) 

J  ^ 


where 


<k  (r)  =  ^  +  i  +!  =  wQ(r)  -  i  (1-x)  r 


(B.3) 


The  boundary  condition  at  x  =  A  is  satisfied  if  Fr(t)  and  Fj  ( r)  are 
identically  zero.  Since  they  are  not,  we  seek  2  NJ  constants,  the  eRj's  and 
e j j ' s  that  minimize  the  error  defined  as 


1 

g  (ej)  r  ^  (Fr2  +  Fj2)  dr. 
0 


(B.4) 


The  minimizing  conditions  are 


3g/9eRj.  =  0, 


ag/aejj  =  0 


j  =  1 , . . . ,  N  J  . 


(B.5) 


Substituting  the  required  terms  into  Eq.  (B.5)  and  performing  the 
manipulations  lead  to  the  linear  system  for  the  e^'s  (j  =  1....NJ): 


f 1  NJ  f 1 

<«Rj  %  *  SIj  *I>  dr  ‘  .  E  .  eRi  t  \  <*Rj  “Ri  +  SIj  “li>  dr]  + 
Jn  1  "  1  Jn 


(B.6a) 


NJ 

E  eIi  C  \  (-w„.  wTl.  +  wT.  w0J  dr] 


i  =  1 


Rj  "Ii  "Ij  "Ri 
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( B .6b ) 


(5Rj  n  ■  5ij  *> dr  - 


<-5Ij  SRi  *  5Rj  *11  >  dr]  * 


NJ  f 1 

.  E  ,  en  [\  (“ij  5n  +  wRj  wri>  dr]- 

1  =  1  Jo 


dj's  are  then  found  from  Eq.  (B.l). 
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LIST  OF  SYMBOLS 


A 

bk 

B 

c 

Cr 

C(LM)P 

C(LIM)P 

C(LSM)P 

C(LSM)PM 

dj 

ej 

f 

F 

h 

h 

j 

k 

KF 

Ko 

K1 

M 


cross-sectional  radius  of  cylinder.  Figure  1  [cm] 

=  c/a,  aspect  ratio  of  cylinder 

bi orthogonal  coefficients  for  series  of  axial 
eigenfunctions,  Eqs.  (2.28)  and  (2.29) 

function  of  r  defined  in  Eq.  (2.26) 

half-height  of  cylinder  [cm] 

natural  oscillation  frequency  of  rotating  liquid/^ 

complex  liquid  pressure  moment  coefficient  =  + 

1  C(LIM)P»  ^1*  (3*3) 

pressure  in-plane  moment  coefficient,  Eq.  (3.18) 
pressure  side  moment  coefficient,  Eq.  (3.18) 
maximum  value  of  C(l$m)p 

coefficient  in  series  given  by  the  3rd  terms  in  Eq. 
(2.17);  see  also  Eq.  (2.34) 

coefficient  defined  in  Eq.  (B.l) 

=  (l-ie)x,  complex  representation  of  angular  motion, 
Eq.  (1.4) 

endwall  boundary  condition  function  of  r,  Eq.  (2.31) 

integral  occurring  in  the  endwall  moment  calculation, 
Eqs.  (3.12)  and  (3.14) 

integral  occurring  in  the  sidewall  moment  calculation, 
Eqs.  (3.8)  and  (3.10) 

index  in  series  given  by  the  3rd  terms  in  Eq.  (2.17) 

index  of  axial  eigenfunction  and  eigenvalue, 

Eqs.  (2.18)  and  (2.19) 

index  of  final  term  of  sin  X^x  series  in  Eq.  (2.17a) 
yaw  amplitude  at  time  t  =  0 

£  T 

=  Kq  e'  9  ,  yaw  amplitude  at  time  t,  Eq.  (1.4) 

=  x  -  V/r 
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LIST  OF  SYMBOLS  (Continued) 


Mj,  M2 

moment  terms,  Eqs.  (3.16)  and  (3.17) 

M(LY)P*  M(LZ)P 

y  and  z  components,  respectively,  of  liquid  moment 
[g  cm2/s2] 

M(LZB)P*  M(LZT)P 

bottom  and  top  wall  contributions,  respectively,  to 
M(LZ)P  cm2/s2] 

M(LZE)P,  M(LZL)P 

endwall  and  sidewall  contributions,  respectively,  to 

m(lz)p  ^  crn2/s2l 

n 

index  of  radial  mode  for  eigenfrequency, 

nYE»  nZE 

components  in  the  y,  z  plane  of  a  unit  vector 
lying  on  the  x-axis 

NO 

number  of  terms  in  the  series  given  by  the  3rd  terms 
in  Eq.  (2.17) 

P 

* 

P 

★ 

pc 

pressure/(pa2$2) 

perturbation  pressure/(KQpa2|2) ,  Eqs.  (2.1)  and  (2.4) 
complex  perturbation  pressure/(Kopa2$2) ,  Eq.  (2.8) 

B  =  Br  +  1  Ei 

r,x  variation  of  perturbation  pressure/(K  pa2|2), 

Eq.  ( 2 . 8d )  0 

Pj(r) 

jth  eigenfunction  in  series  given  by  the  3rd  term  in 

Eq.  ( 2 . 1 7d ) 

P>> 

P>) 

first  term  in  solution  to  p,  Eq.  ( 2 . 17d ) 

coefficient  of  kth  axial  eigenfunction  in  series  given 
by  the  2nd  term  in  Eq.  (2. 17d ) 

P 

★ 

Wedemeyer  model  approximation  to  P  in  core  flow, 

Eqs.  (2.2b)  and  (2.2e) 

★ 

P 

axisymmetric  unperturbed  pressure/(pa2^2) ,  Eq.  (2. Id) 

pl 

P  (r  =  1)  =  constant 

r»r 

(1/a)  x  radial  coordinate  in  inertial  coordinate  system 
and  non-rotating  aeroballistic  system,  respectively 
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LIST  OF  SYMBOLS  (Continued) 


Re 

Reynolds  number  =  a2$/v 

Re£ 

"effective"  Reynolds  number,  Eq.  (2.14) 

t 

time  [s] 

tM 

time  at  which  C(lSM)PM  occurs  for  given  t 

£ 

4>t,  non-dimensional  time,  or  total  angle  of  rotation 

u,u 

(l/[a|3)  x  radial  velocity  component  in  inertial  frame 
and  non-rotating  aeroballistic  frame,  respectively 

v,v 

( 1/ [a $] )  x  azimuthal  velocity  component  in  inertial 

frame  and  non-rotating  aeroballistic  frame, 
respectively 

w,w 

(l/[a$])  x  axial  velocity  component  in  inertial  frame 
and  non-rotating  aeroballistic  frame,  respectively 

★  ★  ★ 

U,  V,  w 

(l/[KQa$])  x  radial,  azimuthal,  axial  perturbation 

velocity  components  in  inertial  system,  Eqs.  (2.1) 
and  (2.4) 

★  *  ★ 

V  V  wc 

(l/[KQa<£])  x  complex  radial,  azimuthal,  and  axial 
perturbation  velocity  components,  Eq.  (2.8) 

U,  V,  w 

(l/[KQa$])  x  r,x  variation  of  perturbation  velocity 
components,  Eqs.  (2.8)  and  (2.9) 

AAA 

V  V  wo 

1st  terms  in  right-hand  sides  of  Eqs.  (2.17a),  (2.17b), 
and  (2.17c) 

v  v  \ 

coefficients  of  sin  xkx  and  cos  X^x  in  series  expansions 
for  u,  y,  w,  g  (2nd  terms  in  Eq.  (2.17)) 

V  V 

coefficients  of  sin  p^x  and  cos  ujX  in  series  expansions 
for  u,  y,  w,  p  (3rd  terms  in  Eq.  (2.17)) 

u,  V,  W 

*  ★  * 

U,  V,  w 

★  *  ★ 

Wedemeyer  core  flow  model  approximations  to  U,  V,  W 

( 1/Ca )  x  radial,  azimuthal,  and  axial  velocity 
components  of  unperturbed  flow,  Eq.  (2.1) 
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LIST  OF  SYMBOLS  (Continued) 


x,  y,  z 

x,  y,  z 


yl»****y6 


rectangular  coordinates  in  inertial  system  (x-axis  along 
trajectory)  [length/a] 

rectangular  coordinates  in  aeroballistic  system 
(x-axis  along  cylinder  axis)  [length/a] 

functions  describing  radial  variation  of  perturbation 
flow  variables,  Eqs.  (2.24)  and  (2.25) 


6c 


6c 


E 


e 


v 

l 

P 

T 


TM 

* 

i>  (r) 


correction  term  in  endwall  boundary  condition  for 
solid  body  rotation 

correction  term  in  ad  hoc  endwall  boundary  condition, 
Eqs.  (2.12)  and  (2.13) 

=  (1/t)  x  yaw  growth  per  radian  of  nutation 

=  0  for  X^  *  0,  =  1  for  x^  =  0 

polar  angles  (azimuthal  coordinates)  in  inertial  and 
aeroballistic  systems,  respectively 

eigenvalue  in  the  axial  problem,  Eqs.  (2.18)  and  (2.19) 

eigenvalues  occurring  in  series  given  by  the  3rd  terms 
in  Eq.  (2.17);  see  discussion  preceding  Eq.  (2.20) 

kinematic  viscosity  of  liquid  [cm2/s] 

vector  describing  angular  motion  of  cylinder,  Eqs.  (1.2) 
and  (1.3) 

density  of  liquid  [g/cm3] 

nutational  frequency  of  cylinder/^ 

t  for  which  C^SMJPM  occurs  at  a  given  t 

spin  rate  of  cylinder  [rad/s],  taken  to  be  positive 

-  wQ  (r)  -  1  (1  -  x)  r,  Eq.  (2.32) 
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ATTN:  B.H.  Little,  Jr. 

G. A.  Pounds 

Dept  72074,  Zone  403 
86  South  Cobb  Drive 
Marietta,  GA  30063 

1  Lockheed  Missiles  and  Space 

Company 

ATTN:  Tech  Info  Center 
3251  Hanover  Street 
Palo  Alto,  CA  94304 

3  Martin-Marietta  Corporation 
ATTN:  S.H.  Maslen 

S.C.  Traugott 

H.  Obremski 
1450  S.  Rolling  Road 
Baltimore,  MO  21227 

2  McDonnell  Douglas  Astronautics 

Corporation 
ATTN:  J.  Xerikos 
H.  Tang 

5301  Bolsa  Avenue 
Huntington  Beach,  CA  92647 


1  Douglas  Aircraft  Company,  Inc. 
ATTN:  T.  Cebeci 
3855  Lakewood  Boulevard 
Long  Beach,  CA  90846 

3  Rockwell  International 

Science  Center 
ATTN:  Or.  V.  Shankar 
Or.  N.  Malmuth 
Dr.  S.  Chakravarthy 
1049  Camino  Dos  Rios 
Thousand  Oaks,  CA  91360 

4  Director 

Sandia  National  Laboratory 
ATTN:  H.W.  Vaughn 

F.G.  Blottner 
W.L.  Oberkampf 
Tech  Library 
Albuquerque,  NM  87115 


No.  of 

Copies  Organization 

2  United  Aircraft  Corporation 
Research  Laboratory 

ATTN:  M.J.  Werle;  Library 

400  Main  Street 

East  Hartford,  CT  06108 

1  LTV  Aerospace  Corporation 
Vought  Systems  Division 
ATTN:  J.M.  Cooksey,  Chief, 

Gas  Dynamics  Lab, 
2-53700 
P.0.  Box  5907 
Dallas,  TX  75222 

1  Arizona  State  University 
Department  of  Mechanical  and 
Energy  Systems  Engineering 
ATTN:  G.P.  Neitzel 
Tempe,  AZ  85281 

1  Cornell  University 
Graduate  School  of 
Aero  Engr 
ATTN:  Library 
Ithaca,  NY  14850 

3  California  Institute  of 

Technology 
ATTN:  Tech  Library 

H.B.  Keller,  Math  Dept 
D.  Coles,  Aero  Dept 
Pasadena,  CA  91109 

1  Illinois  Institute  of  Technology 
ATTN:  H.  M.  Nagib 
3300  South  Federal 
Chicago,  IL  60616 

1  The  Johns  Hopkins  University 
Dept  of  Mech  and  Materials  Sci. 
ATTN:  S.  Corrsin 
Baltimore,  MD  21218 
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4  Director 

Johns  Hopkins  University 
Applied  Physics  Laboratory 
ATTN:  Dr.  R.D.  Whiting 
Dr.  D.A.  Hurd if 

Dr.  R.S.  Hirsh 
Mr.  E.R.  Bohn 
Johns  Hopkins  Road 
Laurel,  MD  20707 

3  Massachusetts  Institute  of 
Technology 
ATTN:  E.  Covert 

H.  Greenspan 
Tech  Lib 

77  Massachusetts  Avenue 
Cambridge,  MA  02139 

2  North  Carolina  State  Univ 
Mechanical  and  Aerospace 
Engineering  Department 
ATTN:  F.F.  DeJarnette 
J.C.  Williams 
Raleigh,  NC  27607 

1  Northwestern  University 
Department  of  Engineering 
Science  and  Applied 
Mathematics 
ATTN:  Dr.  S.H.  Davis 
Evanston,  IL  60201 

1  Notre  Dame  University 
Department  of  Aero  Engr 
ATTN:  T.J.  Mueller 
Notre  Dame,  IN  46556 

2  Ohio  State  University 
Dept  of  Aeronautical  and 

Astronautical  Engineering 
ATTN:  S.L.  Petrie 

O.R.  Burggraf 
Columbus,  OH  43210 


No.  of 

Copies  Organization 

2  Polytechnic  Institute  of 

New  York 

ATTN:  G.  Moretti 
S.G.  Rubin 
Route  110 

Farmingdale,  NY  11735 

3  Director 

Forrestal  Research  Ctr 
Gas  Dynamics  Laboratory 
Princeton  University 
ATTN:  S.M.  Bogdonoff;  S.I.  Cheng 
Tech  Library 
Princeton,  NJ  08540 

1  Purdue  University 

Thermal  Science  &  Prop  Ctr 
ATTN:  Tech  Library 
W.  Lafayette,  IN  47906 

1  Rensselaer  Polytechnic 
Institute 

Department  of  Math  Sciences 
ATTN:  R.C.  Di  Prima 
Troy,  NY  12181 

1  Southern  Methodist  University 
Dept,  of  Civil  and  Mechanical 
Engineering 
ATTN:  R.  L.  Simpson 
Dallas,  TX  75222 


1  San  Diego  State  University 
Department  of  Aerospace  Engr 
and  Engineering  Mechanics 
College  of  Engineering 
ATTN:  K.C.  Wang 
San  Diego,  CA  92115 
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1  Southwest  Research  Institute 
Applied  Mechanics  Reviews 
8500  Culebra  Road 

San  Antonio,  TX  78228 

2  Stanford  University 

Dept  of  Aeronautics/Astronautics 
ATTN:  Dr.  J.L.  Steger 

Dr.  S.  Chakravarthy 
Stanford,  CA  94305 

1  Texas  A&M  University 
College  of  Engineering 
ATTN:  R.H.  Page 
College  Station,  TX  77843 

1  University  of  California  - 
Davis 

ATTN:  H.A.  Dwyer 
Davis,  CA  95616 

1  University  of  California  - 

Berkeley 

Department  of  Aerospace 
Engineering 
ATTN:  M.  Holt 
Berkeley,  CA  94720 

2  University  of  California  - 

San  Diego 

Department  of  Aerospace 
Engineering  and  Mechanical 
Engineering  Sciences 
ATTN:  P.  Libby 

Tech  Library 
La  Jolla,  CA  92093 

1  University  of  Cincinnati 
Department  of  Aerospace 
Engi neeri ng 
ATTN:  R.T.  Davis 
Cincinnati,  OH  45221 

1  University  of  Colorado 

Department  of  Astro-Geophysics 
ATTN:  E.R.  Benton 
Boulder,  CO  80302 


2  University  of  Maryland 
ATTN:  W.  Melnik 

J.D.  Anderson 
College  Park,  MD  20740 

1  University  of  Maryland  - 
Baltimore  County 
Department  of  Mathematics 
ATTN:  Dr.  Y.M.  Lynn 
5401  Wilkens  Avenue 
Baltimore,  MD  21228 

1  University  of  Santa  Clara 
Department  of  Physics 
ATTN:  R.  Greeley 

Santa  Clara,  CA  95053 

2  University  of  Southern 

California 

Department  of  Aerospace 
Engineering 
ATTN:  T.  Maxworthy 
P.  Weidman 

Los  Angeles,  CA  90007 

2  University  of  Michigan 
Department  of  Aeronautical 
Engi neeri ng 
ATTN:  W.W.  Wilmarth 
Tech  Library 

East  Engineering  Building 
Ann  Arbor,  MI  48104 

2  University  of  Rochester 
Department  of  Mechanical 
and  Aerospace  Sciences 
ATTN:  R.  Gans 

A.  Clark,  Jr. 
Rochester,  NY  14627 

1  University  of  Tennessee 
Department  of  Physics 
ATTN:  Prof.  W.E.  Scott 
Knoxville,  TN  37916 
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1  University  of  Texas 
Department  of  Aerospace 
Engineering 

ATTN:  J.C.  Westkaemper 
Austin,  TX  78712 

1  University  of  Virginia 
Department  of  Aerospace 
Engineering  &  Engineering 
Physics 

ATTN:  I.D.  Jacobson 
Charlottesville,  VA  22904 

1  University  of  Virginia  Research 
Laboratories  for  the 
Engineering  Sciences 
ATTN:  Prof.  H.  G.  Wood 
P.0.  Box  3366 
University  Station 
Charlottesville,  VA  22903 

1  University  of  Washington 
Department  of  Mechanical 
Engineering 
ATTN:  Tech  Library 
Seattle,  WA  98105 

1  University  of  Wyoming 
ATTN:  D.L.  Boyer 
University  Station 
Laramie,  WY  82071 

1  U.S.  Military  Academy 
Department  of  Physics 
ATTN:  MAJ  G.  Heuser 
West  Point,  NY  10996 


No.  of' 
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1  Woods  Hole  Oceanographic 
Institute 

ATTN:  J.A.  Whitehead 
Woods  Hole,  MA  02543 

3  Virginia  Polytechnic  Institute 
and  State  University 
Department  of  Aerospace 
Engineering 
ATTN:  Tech  Library 
Dr.  W.  Saric 
Dr.  T.  Herbert 
Blacksburg,  VA  24061 


Aberdeen  Proving  Ground 

Director,  USAMSAA 
ATTN:  DRXSY-D 

DRXSY-MP,  H.  Cohen 

Commander,  USATECOM 
ATTN:  DRSTE-TO-F 


Commander,  CRDC,  AMCCOM 
ATTN:  DRSMC-CLN 
W.  C.  Dee 
DRSMC-CLB-PA 
M.  C.  Miller 
DRSMC-CLJ-L 
DRSMC-CLB-PA 
DRSMC-CLN 
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USER  EVALUATION  OF  REPORT 


Please  take  a  few  minutes  to  answer  the  questions  below;  tear  out 
this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  place 
in  the  mail.  Your  comments  will  provide  us  with  information  for 
improving  future  reports. 

1 .  BRL  Report  Number _ _ 

2.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related 
project,  or  other  area  of  interest  for  which  report  will  be  used.) 


3.  How,  specifically,  is  the  report  being  used?  (Information 
source,  design  data  or  procedure,  management  procedure,  source  of 
ideas,  etc.) _ 


4.  Has  the  information  in  this  report  led  to  any  quantitative 
savings  as  far  as  man-hours/contract  dollars  saved,  operating  costs 
avoided,  efficiencies  achieved,  etc.?  If  so,  please  elaborate. 


5.  General  Comments  (Indicate  what  you  think  should  be  changed  to 
make  this  report  and  future  reports  of  this  type  more  responsive 
to  your  needs,  more  usable,  improve  readability,  etc.) _ 


6.  If  you  would  like  to  be  contacted  by  the  personnel  who  prepared 
this  report  to  raise  specific  questions  or  discuss  the  topic, 
please  fill  in  the  following  information. 

Name : 


Telephone  Number: 
Organization  Address: 


I 

i 


L 


